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Complex  Variable  Methods  for  Fatigue  Sensitivity  Analysis 

A.  Voorhees1,  R.  Bagley2,  H.  Millwater3, 

Graduate  student,  2Professor,  3Associate  Professor 
Department  of  Mechanical  Engineering,  Univ.  of  Texas  at  San  Antonio 

Abstract 

Complex  variable,  numerical  differentiation  techniques  have  proven  useful  in  many  fields  of  engineering  analysis. 
Complex  Taylor  series  expansion  and  Fourier  differentiation  are  two  such  complex  variable  methods.  This  paper  adapts 
the  use  of  both  complex  Taylor  series  expansion  and  Fourier  differentiation  for  fatigue  sensitivity  analysis.  The 
sensitivity  of  the  number  of  cycles  to  failure  to  input  parameters  and  initial  conditions  has  been  determined  by  traditional 
numerical  differentiation  techniques  as  well  as  through  complex  variable  methods.  Both  complex  Taylor  series  expansion 
and  Fourier  differentiation  have  been  found  to  be  more  accurate  and  stable  than  traditional  central  differencing. 

Introduction 

The  fatigue  of  aircraft  components  due  to  repeated  loading  cycles  is  a  costly  problem  that  may  also  poses  safety 
risks  to  pilots  and  crewmembers.  If  a  part  fails  due  to  fracture  it  can  mean  costly  repairs  or  even  retirement  of  the  plane. 
The  design  and  maintenance  of  aircraft  should  take  into  account  the  effects  of  fatigue.  The  goal  of  proper  design  and 
maintenance  should  be  to  increase  the  life  of  the  aircraft  while  ensuring  its  safety.  At  the  same  time  it  is  desirable  to 
reduce  the  cost  of  operation  and  maintenance.  In  order  to  determine  the  optimal  design  and  maintenance  procedures,  it  is 
necessary  to  know  the  relative  importance  of  the  variables  in  the  problem.  This  relative  importance  can  be  measured  by 
the  sensitivity  or  the  partial  derivative.  If  the  sensitivity  information  was  known,  it  would  become  possible  to  run  an 
optimization  routine  to  determine  the  optimal  dimensions,  or  material  properties  of  a  part.  It  would  also  be  possible  to 
determine  the  optimal  parameters  of  an  inspection  routine  designed  to  identify  small  cracks  before  they  threaten  the  life  of 
the  part.  Traditional  numerical  differentiation  should  allow  for  the  calculation  of  such  sensitivities,  however  newer 
complex  variable  methods  promise  more  accuracy,  especially  in  the  calculation  of  higher  order  derivatives. 

Numerical  differentiation  techniques  such  as  finite  differencing  are  commonly  used  to  evaluate  the  derivatives  of 
implicit  functions.  In  the  simplest  form  of  finite  differencing  the  derivative  of  a  function  at  a  point  can  be  calculated  by 
evaluating  the  function  at  a  second  nearby  point  located  along  the  real  axis.  This  nearby  point  is  the  sample  point,  and  the 
distance  between  the  sample  point  and  the  original  point  is  the  sampling  radius.  The  difference  in  the  value  of  the 
function  evaluated  at  the  sample  point  and  at  the  original  point  divided  by  the  sampling  radius  yields  an  estimate  of  the 
first  derivative.  Higher  order  derivatives  as  well  as  more  theoretically  accurate  derivatives  can  be  calculated  by 
evaluating  the  function  at  more  sample  points  and  performing  more  difference  operations.  The  finite  difference  method  is 
easy  to  perform,  however  it  is  only  accurate  when  the  difference  between  the  sample  points  is  small.  When  subtracting 
small  numbers,  machine  round  off  can  also  introduce  error  meaning  that  there  is  a  limit  to  the  accuracy  of  the  method. 
This  becomes  more  of  a  problem  when  calculating  higher  order  derivatives,  since  more  subtraction  operations  are 
required. 

An  alternative  to  finite  differencing  is  complex  Taylor  series  expansion  (CTSE).  CTSE  was  originally  described  by 
Lyness  and  Moler  in  the  1960’s  but  wasn’t  brought  to  the  attention  of  the  engineering  community  until  the  1990’s  when 
Squire  showed  the  accuracy  of  the  method  for  computing  the  first  derivative  of  several  explicit  functions  [1]  [2]  [3] .  In 
CTSE,  the  sample  point  is  complex.  The  first  derivative  is  calculated  simply  by  determining  the  imaginary  component  of 
the  function  evaluated  at  the  sample  point.  For  the  first  derivative  CTSE  involves  no  difference  operation  and  the 
theoretical  error  is  smaller  than  for  finite  differencing.  In  order  to  calculate  the  higher  order  derivatives,  additional  sample 
points  need  to  be  taken  and  more  difference  operations  become  necessary. 

Fourier  Differentiation  (FD)  is  another  alternative  method  using  sampling  in  the  complex  plane.  Fourier 
differentiation  requires  the  evaluation  of  several  sample  points  along  a  circular  contour  around  the  initial  point  in  the 
complex  plane.  This  technique  was  first  described  by  Lyness  et  al  in  1967  [1][2].  The  function  is  evaluated  at  these 
sample  points  and  a  fast  Fourier  transform  (FFT)  of  the  evaluated  sample  points  is  used  to  calculate  the  derivative.  The 
use  of  the  FFT  as  a  method  to  calculate  the  derivatives  was  described  by  Lyness  et  al  [4] .  More  recently,  FD  was  shown 
to  be  useful  in  the  computation  of  senstivities  of  implicit  functions  [5].  One  benefit  of  FD  is  that  the  results  of  the  FFT 
also  contain  higher  order  derivatives.  The  number  of  derivatives  that  can  be  obtained  from  the  FFT  is  related  to  the 
number  of  sample  points  chosen.  Furthermore,  FD  is  theoretically  more  accurate  than  CTSE  and  finite  differencing. 
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CTSE  has  already  caught  on  in  many  engineering  fields.  In  fluid  dynamics,  CTSE  has  been  used  to  find  sensitivities 
for  the  solution  of  the  Navier-Stokes  equation  [6].  Furthermore  researchers  have  been  able  to  apply  CTSE  techniques  to 
finite  element  methods  in  the  field  of  aerodynamics  and  aero -structural  analysis  [7]  [8].  CTSE  has  also  proven  to  be  a 
valuable  technique  in  the  study  of  heat  transfer  [9],  dynamic  system  optimization  [10],  and  has  even  been  incorporated 
into  pseudospectral  [11]  and  eigenvalue  sensitivity  methods  [12].  A  thorough  review  of  the  literature  has  not  turned  up 
any  instances  of  CTSE  being  applied  to  fatigue  problems.  Standard  finite  difference  techniques  have  been  used  in  fatigue 
problems  [13]  [14],  as  well  as  probabilistic  sensitivity  methods  [15].  By  adopting  complex  variable  numerical 
differentiation  techniques,  it  should  be  possible  to  obtain  more  accurate  sensitivity  estimates. 

The  goal  of  this  paper  is  to  explore  the  use  of  complex  variable  methods  to  determine  sensitivities  of  the  cycles -to- 
failure  with  respect  to  input  quantities  such  as  loading  amplitude,  crack  growth  parameters,  initial  crack  size,  fracture 
toughness,  etc. 

Additionally  higher  order  derivatives  will  be  calculated  which  may  be  used  to  create  Taylor  series  estimates  of  the 
number  of  cycles  to  failure  as  a  function  of  the  aforementioned  initial  conditions  and  properties. 


II.  Methodology 

A.  Numerical  Differentiation 

The  derivative  of  a  real  function  can  be  defined  as  the  instantaneous  change  in  the  value  of  the  function  at  a  given 
location.  The  limit  definition  of  the  derivative  states: 


/'(*„)  =  lim; 


(1) 


x-x0 

where  f’(x0)  is  the  derivative  of  f(x0).  The  use  of  this  definition  requires  f(x)  to  be  continuous  at  x0.  In  situations  where 
this  limit  can  not  be  evaluated  easily  or  exactly  (i.e.  when  f(x)  is  implicit)  f’(x0)  can  be  calculated  by  numerical 
approximation  of  the  limit.  The  most  traditional  approximation  method  is  finite  differencing.  The  finite  differencing 
method  approximates  the  derivative  by  setting  f(x)  in  equation  1  to  f(xQ+h)  where  h  is  very  small. 
f(x0  +  h)-f(x0) 


/'(*„) 


h 


(2) 


This  approach  requires  the  evaluation  of  the  function  at  two  points,  the  initial  point  and  the  sample  point  x0+h.  When 
h  is  positive,  the  method  is  referred  to  as  forward  differencing  and  when  h  is  negative  it  is  referred  to  as  backwards 
differencing.  A  more  accurate  approximation  can  be  obtained  by  averaging  the  forward  and  backwards  differences  for 
equal  values  of  h. 

fjx0)-f{x0-h)  |  f(xo  +  h)-f(xo)  f(x0  +  h)-f(x0-h) 

2  h  2  h  2  h 


(3) 


This  averaging  technique  is  called  central  differencing,  and  is  often  the  method  of  choice  in  numerical  differentiation. 
The  increased  accuracy  of  central  differencing  arises  from  the  cancellation  of  even,  higher  ordered  terms  in  the  Taylor 
series.  This  pushes  the  accuracy  from  0(h)  for  forward  differencing  to  0(h2)  for  central.  In  this  paper,  central 
differencing  was  selected  for  comparison  with  CTSE  and  FD.  Figure  l.A  displays  the  location  of  sampling  points  in  the 
complex  plane. 

Finite  differencing  can  also  be  used  to  evaluate  higher  order  derivatives.  The  second  derivative  of  f(x)  can  be 
defined  as  the  derivative  of  a  derivative. 


f2\x0)=  lim 


lim 


/(*2)  "A*.)  ,;/(*,)-/«> 


lim 

X— 


(4) 


xl~Xo 


Approximating  the  limits  by  sampling /(A)  at  f(xa),f(x0+h),  and  f(x0+2h)  where  h  is  small  the  second  forward  derivative  is. 
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Figure  1.  The  Location  of  Sampling  Point  in  the  Complex  Plane.  A.  The  location  of  the  first  4 
sampling  points  centered  around  the  initial  point  in  central  differencing.  B.  The  location  of  the  first 
two  sample  points  in  CTSE.  C.  8  sample  points  for  FD. 


f(2\xo) 


f(x0  +  2h)-f(x0  +  h)  f(x0  +  h)-f(x0) 


(5) 


f(2\x0) 


f(xe  +  2h)-2f(xe  +  h)  +  f(xe) 

h 2 


Note  here  that  two  sample  points  are  needed  to  calculate  the  second  derivative.  The  central  difference  for  the  second 
derivative  is  calculated  by  sampling /(v)  at  f(x0-h),f(x0),  and  f(x0+h). 

/( 2)(x  )  /(*„  +  h)~  2/(x0)  +  f(x0  -  h)  (6) 

h2 

Finite  differencing  techniques  can  be  extended  to  higher  order  derivatives  still.  Unfortunately  this  also  requires 
additional  sample  points.  When  the  number  of  sample  points  increases  the  error  can  become  quite  large.  This  error  comes 
from  two  different  sources.  The  first  source  is  subtractive  cancellations  in  which  the  difference  between  two  sample 
points  becomes  too  small  to  be  accurately  calculated.  This  type  of  error  occurs  mainly  in  higher  order  derivatives  when 
the  step  size  is  too  small.  The  second  source  of  error  comes  from  the  approximation  itself.  This  places  the  user  at  a  crux, 
as  h  must  be  selected  so  as  to  be  neither  too  small  nor  too  big. 

By  sampling  in  the  complex  plane  rather  than  solely  along  the  real  axis,  some  of  these  problems  can  be  remedied. 
Derivatives  can  also  be  defined  in  the  complex  plane  using  limits 

f\za)  =  lim  (7) 

z_>z0  z  -  z„ 
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where  z  is  the  complex  variable.  In  order  for  a  function  to  have  a  derivative  at  z  =  z0  it  must  be  analytic  in  the 
neighborhood  surrounding  z0-  A  function  is  analytic  if  its  derivative  is  identical  no  matter  what  direction  in  the  complex 
plane  z  approaches  from. 

If  the  complex  function  is  defined  as 

f(z)  =  f(x  +  i-y)  =  u(x,y )  +  i  ■  v(x,y) ,  (8) 

then  the  derivative  can  be  written  as 


f\z)  = 


di(x,y)  .  dv(x,y)  di(x,y)  .  dv(x,y) 


-  +  l 


dy 


+  l  ' 


dy 


(9) 


dx  dx 

The  Cauchy-Riemann  equations  state  that 

da  _  dv  du  _  dv 

dx  dy  dy  dx ’ 

which  allows  equation  9  to  be  rewritten  as 

....  da  .di  dv  ,d> 

f  (z)  =  —  +  i —  = - 1  — . 

dx  dy  dy  dx 

When  the  function  of  interest  is  real  valued,  only  the  real  part  of  f’(z)  need  be  considered. 

real(f'(z))  =  f(x)  =  ^-=^. 

ax  ay 

This  reveals  that  the  real  part  of  the  first  derivative  can  be  determined  by  examining  the  change  along  either  the  real  axis 
(v)  or  the  imaginary  axis  (y).  The  limit  defined  by  the  partial  of  u  with  respect  to  x  is  the  same  as  the  one  approximated 
by  finite  differencing. 


(10) 


(11) 


(12) 


P(x  )  _  _  lim  u(x,y0)  -  u(x0,y0)  v  -  a(*0,T0) 

dx  x-x„  x-x„ 


(13) 


Approximating  the  limit  defined  by  the  partial  of  v  with  respect  to  y  leads  to  a  whole  new  type  of  numerical 
differentiation. 


f\x  )  =  —  =nm  .. 

^  ^  j-j,,  y-ya 


(14) 


Replacing  y  with  h,  and  accounting  for  the  fact  that  y0  and  v(x0,y0)  are  zero  due  to  the  fact  that  the  initial  point  is  real 
valued,  equation  14  becomes. 

v(xa,h)  _  imag{f(x0  +  ih )) 


f(x0) 


(15) 
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■  + . 


h  h 

In  this  manner,  the  derivative  can  be  approximated  by  adding  a  small  imaginary  component  to  x0.  this  method  is  referred 
to  as  CTSE  and  it  gets  its  name  due  to  the  fact  that  the  method  is  easily  derived  from  the  complex  Taylor  series 

/(*. +m= /w + fH*.)  ■  '-j-  -  -  /<J,o„)  ■,'h> 

Taking  the  imaginary  part  of  the  Taylor  series  results  in 

imag(f(x0  +  ih ))  =  f(l\x0 )  •  ^  -  f°\x0)  ~ ... 

which  is  solved  for  T 1  ‘f  x,,) 


(16) 


(17) 


.... 


(18) 


h  °  3! 

It  can  be  noted  that  the  error  due  to  the  approximation  is  still  0(h2),  however  there  is  no  difference  equation  in  the  first 
derivative,  meaning  that  there  is  no  lower  limit  to  the  selected  step  size. 

CTSE  can  also  be  used  to  determine  higher  order  derivatives.  Taking  the  real  part  of  the  Taylor  series  and  rearranging 
will  allow  for  the  calculation  of  the  second  derivative. 
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real(f(x.  +  ft))  =  /(*«,)  -/mW  ■  ^  +  f\xj  ~ ... 


li 2  A4 

"  r(4)  .  « 


2! 


4! 


U 


(19) 


ji 2)(x  ^  2(/(x0)-rgfl/(/(x0 +  //?)))  2  ^(4)^  ^ 


62  4! 

In  order  to  obtain  estimates  for  the  higher  derivatives,  additional  sampling  points  must  be  taken.  The  Taylor  series  of 
each  sample  point  can  be  coupled  to  create  a  system  of  equations,  that  can  be  solved  for  the  desired  derivatives  as  seen  in 
equation  20. 

imag(f(x0  +  ih))  =  /'>(*„)  •  £  -  fm(xj  ■  ^  +  f!\x„) 

1  3  3  3  5  5  (20) 
imag(f(x0  +  i2h))  =  /">(*,)  ■  “  -  /">(*„)  •  ?2L  +  /»(*„) .  ... 

It  is  clearly  seen  from  equations  19  and  20  that  the  calculation  of  the  second  order  and  higher  derivatives  requires 
difference  operations,  and  as  such  the  method  will  face  problems  similar  to  those  of  finite  differencing  in  that  h  must  not 
be  too  small  nor  too  big.  CTSE  does  have  the  advantage  that  only  one  sample  point  is  needed  to  calculate  both  the  first 
and  second  derivatives,  however.  Central  differencing  requires  two  sample  points  to  obtain  the  first  two  derivatives. 
Sampling  need  not  be  limited  solely  to  the  axes  though.  Sample  points  can  be  selected  at  any  point  in  the  complex  plane. 
For  instance  in  FD  sample  points  are  located  along  a  contour  encircling  the  initial  point  as  seen  in  figure  l.C.  FD  is  a 
numerical  technique  based  on  the  Cauchy  integral  theorem. 

/(”)(z)  =  — j  ^  ,  4g  (21) 

2m^  (<f-z)  1 

Derivatives  of  any  order  can  thus  be  calculated  by  evaluating  a  simple  contour  integral.  One  way  to  evaluate  this  integral 
is  to  take  advantage  of  its  similarities  to  the  Fourier  integral.  Through  a  change  in  variables, 

z  — »  zo,  and  ^  —  zo  +  cew  (22) 

the  Cauchy  integral  is  transformed  into  a  Fourier  integral 


fn\z)  =  —  \2\f(zo  +  cew)el,,od0. 
Inc11  J° 


(23) 


This  means  that  the  function /(z)  can  be  evaluated  by  using  Fourier  transforms.  If  the  function  in  the  Cauchy  integral 
theorem /(z)  is  analytic  in  a  circular  neighborhood  centered  on  z0  with  a  radius  of  h,  then  its  Taylor  series  exists  and 
converges  within  that  neighborhood.  The  resulting  Taylor  series  is 


fiz^),±rMc^ 

— a  n\ 


n  in6 


(24) 


n=  0 


Notice  that  this  function  is  periodic  in  Shaving  period  2n.  Also  note  that  each  term  in  the  series  oscillates  at  positive 
integer  multiples  of  0.  Because  of  this  pattern  of  oscillation,  this  Taylor  series  is  seen  to  take  on  the  properties  of  a 
periodic  Fourier  series.  In  addition  each  term  in  the  Taylor  series  is  orthogonal  to  all  other  terms  on  the  interval  0  to  2n. 


N 


/(z,W“)  =  2>„e“’ 


(25) 


n=  0 


Sampling  the  function  f(z0+hel6)  at  several  values  of  0  between  0  and  2n,  and  then  taking  the  FFT  of  the  sampled  data, 
allows  for  the  calculation  of  the  Fourier  coefficients  an.  These  coefficients  can  then  be  related  back  to  the  Taylor  series 
terms,  and  hence  the  derivatives  of  the  function. 

f^n\z  )  a 

- - ^  =  —  (26) 

n\  c 

The  number  of  derivatives  that  can  be  calculated  is  directly  related  to  the  number  of  sample  points  chosen.  If  N 
sample  points  are  selected,  then  the  first  N/2  terms  in  the  Taylor  series  can  be  calculated.  Due  to  the  orthogonality  of  the 
Fourier  series  terms,  the  order  the  estimation  error  is  also  related  to  the  number  of  sample  points  chosen.  For  example, 
when  the  first  8  terms  of  the  Taylor  series  are  computed,  then  the  estimation  error  will  be  0(h8)  for  every  derivative.  This 
allows  for  the  use  of  much  bigger  sampling  radii,  at  least  an  order  of  magnitude  larger,  than  the  one’s  employed  in  central 
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differencing  and  CTSE.  Unfortunately,  error  occurs  in  the  FFT  when  the  sampling  radius  becomes  too  small,  just  as  with 
the  other  two  methods. 

B.  Fatigue  Analysis 

A  bilinear  Paris  Faw  relationship  was  chosen  to  model  the  da/dN  vs.  AK  curve.  Material  properties  for  the  model 
were  taken  from  experimental  data  for  titanium  6-4.  The  crack  was  chosen  to  be  a  small,  two-dimensional,  centered 
surface  crack.  Variable  amplitude  loading  was  included  in  the  model,  and  again  experimental  data  was  used  for  the  load 
history.  A  forward  Eulerian  differential  equation  solver  was  written  in  Matlab  to  solve  the  crack  growth  equation,  as  this 
is  the  standard  solution  method  in  the  fatigue  field.  This  differential  equation  solver  returned  the  number  of  cycles  to 
failure  as  well  as  the  final  length  and  depth  of  the  crack. 

III.  Results 

It  is  assumed  in  numerical  differentiation  that  the  function  being  differentiated  and  its  derivatives  are  smooth 
functions.  This  is  likely  not  the  case  for  fatigue  analysis  since  the  function  being  differentiated  is  a  first  order  numerical 
solution  to  a  differential  equation.  This  lack  of  smoothness  will  prevent  the  accurate  calculation  of  derivatives  of  an 
arbitrary  order.  It  may  be  possible  to  calculate  derivatives  up  to  a  certain  order  if  the  numerical  solution  and  its  first  few 
derivatives  are  nearly  smooth  over  a  local  region  around  the  initial  point  of  differentiation.  In  order  to  test  the  validity  and 
performance  of  the  three  numerical  differentiation  techniques,  first  and  second  order  derivatives  were  calculated  at  several 
points  in  a  small  neighborhood.  The  differential  equation  solver  allows  the  user  to  specify  the  maximum  percent  growth 
in  crack  length  or  depth  per  solution  step,  referred  to  as  the  step  size  of  the  differential  equation  solver.  Decreasing  the 
step  size  increases  the  amount  of  time  required  to  complete  the  analysis  but  it  results  in  smoother  and  more  accurate 
solutions.  Figure  2  displays  the  first  order  sensitivity  of  the  number  of  cycles  to  failure  with  respect  to  the  initial  crack 
depth  as  calculated  by  CD,  with  h  =  0.0001,  and  figure  3  displays  the  second  order  sensitivity  calculated  by  CD  for  three 
different  step  sizes:  2%,  1%,  and  0.1%. 


Figure  2.  The  First  Derivative  of  the  Number  of  Cycles  to  Failure  with  Respect  to  the  Initial  Crack 
Depth  for  3  Different  Step  Sizes  Calculated  by  CD 
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Second  Derivative:  Step  Size  =  0.02 


Second  Derivative:  Step  Size  =  0.01 


45  46  47  48  49  5  5t  52  53  54  55 


Figure  3.  The  Second  Derivative  of  the  Number  of  Cycles  to  Failure  with  Respect  to  the  Initial  Crack 
Depth  for  3  Different  Step  Sizes  Calculated  by  CD 


From  these  results  it  is  shown  that  first  and  second  order  derivatives  can  be  calculated  smoothly  across  a  local  range  of  the 
input  variable.  It  is  also  shown  that  decreasing  the  step  size  increases  the  smoothness  of  the  second  derivative  and  allows 
for  accurate  calculation  of  the  derivatives.  The  first  and  second  order  sensitivities  calculated  by  CTSE  with  h=0.0001 ,  are 
shown  in  figures  4  and  5  respectively. 

These  results  show  more  stability  at  every  differential  equation  step  size  than  the  derivatives  calculated  via  CD  for 
the  second  derivative,  while  maintaining  stable  first  derivatives.  Derivatives  were  also  calculated  using  a  four  point  FD 
method  with  the  sampling  radius  c=0.0001.  For  the  four  point  FD  method  the  operations  performed  by  the  FFT  result  in 
derivatives  that  are  identical  to  the  average  of  the  CTSE  and  CD  results.  As  such  the  behavior  of  the  sensitivities  over  a 
local  neighborhood  is  similar  to  that  of  CD.  The  results  for  FD  are  not  pictured.  Sensitivities  were  also  calculated  with 
respect  to  several  other  parameters  such  as  material  properties  like  the  Paris  constant,  and  the  maximum  applied  stress.  In 
all  of  these  cases  CTSE  returned  smoother  first  and  second  order  derivatives. 

The  real  advantage  of  FD  over  CTSE  and  CD  is  that  it  allows  for  the  accurate  calculation  of  higher  order  derivatives 
if  they  are  available,  i.e.  smooth.  Using  an  8  point  FD  method  with  c=0.0003 ,  the  first  6  sensitivities  of  the  number  of 
cycles  to  failure  were  calculated  with  respect  to  the  initial  crack  length.  These  results  are  shown  for  a  differential 
equation  step  size  of  1%  in  figure  6  and  .1%  in  figure  7. 
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Figure  4.  The  First  Derivative  of  the  Number  of  Cycles  to  Failure  with  Respect  to  the  Initial  Crack 
Depth  for  3  Different  Step  Sizes  Calculated  by  CTSE 


aO  *  ioJ 

Figure  5.  The  Second  Derivative  of  the  Number  of  Cycles  to  Failure  with  Respect  to  the  Initial  Crack 
Depth  for  3  Different  Step  Sizes  Calculated  by  CTSE 
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1st  Derivative  of  Initial  Crack  Length 


x  l0<<  2nd  Derivative  of  Initial  Crack  Length 


5th  Derivative  of  Initial  Crack  Length 


Figure  6.  The  First  Six  Derivatives  of  the  Number  of  Cycles  to  Failure  With  Respect  to  the  Initial 
Crack  Length  as  Calculated  by  8  point  FD.  For  a  Differential  Equation  Solver  Step  Size  of  1% 


It  is  seen  that  for  differential  equation  solver  step  size  of  1%  that  only  the  first  order  sensitivity  displays  good  stability. 
However  reducing  the  step  size  dramatically  increases  the  smoothness  of  the  higher  order  derivatives.  With  enough 
computational  effort  it  may  be  possible  to  calculate  the  higher  order  derivatives  with  good  stability  over  the  local  domain. 
It  is  also  important  to  note  that  the  sensitivities  with  respect  to  the  initial  crack  length  show  the  least  stability.  Although 
not  shown,  sensitivities  with  respect  to  the  Paris  constants  showed  excellent  stability  through  the  sixth  derivative  at  a 
differential  equation  solver  step  size  of  1%. 
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1st  Derivative  of  Initial  Crack  Length 


2nd  Derivative  of  Initial  Crack  Length 


3rd  Derivative  of  Initial  Crack  Length 
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4th  Derivative  of  Initial  Crack  Length 


x  10-7  6th  Derivative  of  Initial  Crack  Length 
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Figure  7.  The  First  Six  Derivatives  of  the  Number  of  Cycles  to  Failure  With  Respect  to  the  Initial 
Crack  Length  as  Calculated  by  8  point  FD.  For  a  Differential  Equation  Solver  Step  Size  of  0.1% 


IV.  Discussion  and  Conclusions 

Complex  variable  methods  numerical  differentiation  techniques  can  be  used  to  calculate  sensitivities  useful  for 
fatigue  analysis.  CTSE  outperforms  the  traditional  CD  method  for  the  calculation  of  first  and  second  order  sensitivities. 
In  addition  to  being  more  stable  across  a  local  domain  including  the  initial  point,  the  CTSE  method  also  requires  fewer 
sampling  points,  which  may  increase  calculation  speed  depending  on  the  fatigue  analysis  code. 

All  numerical  differentiation  methods  showed  more  stability  in  the  derivatives  when  the  step  size  of  the  differential 
equation  solver  was  decreased.  Decreasing  the  step  size  increases  the  smoothness  of  the  function  and  its  derivatives, 
which  allows  for  the  more  accurate  calculation  of  derivatives.  This  result  also  hints  that  the  use  of  more  accurate 
differential  equation  solver  such  as  a  fourth  order  Runge-Kutta  scheme  may  increase  the  accuracy  and  stability  of  the 
numerical  differentiation  techniques. 

FD  does  not  outperform  CTSE  for  the  calculation  of  low-order  sensitivities  but  it  is  a  useful  tool  in  the  calculation  of 
higher  order  derivatives.  FD  can  be  used  to  calculate  smooth  sixth  order  sensitivities  for  some  input  variables.  For  cases 
where  the  high  order  sensitivities  are  not  very  stable,  decreasing  the  step  size  of  the  differential  equation  solver  increases 
the  stability.  Thus  given  enough  computational  power,  high  order  derivatives  for  any  input  variable  can  be  calculated  via 
FD. 
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